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ABSTRACT 

Recently  acquired  high  pressure  ablation  data  for  graphite  and  carbon 
phenolic  are  analyzed  with  the  objective  of  establishing  the  magnitude  of 
discrepancies  between  measured  ablation  rates  and  sirtooth-wall  "thermochemical 
only"  predictions.  Data  were  obtained  under  the  SAMSO  sponsored  Thermal  Pro¬ 
tection  Program  and  are  believed  of  sufficient  quality  to  reach  meaningful  con¬ 
clusions  with  respect  to  the  magnitude  of  ablation  enhancement  resulting  from 
the  combined  effects  of  rough-wall  and  mechanical  erosion  phenomena.  Compari¬ 
sons  between  predicted  and  measured  recession  are  presented  for  laminar  and 
turbulent  conditions  for  models  tested  in  the  Cornell  Wave  Superheater,  The 
AFFDL  50  MW  arc- jet,  and  the  AEDC  ballistic  range.  The  ratio  of  predicted- to- 
measured  recession  for  laminar  and  turbulent  flow  for  both  materials  correlate 
moderately  well  with  local  pressure. 
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SECiaON  1 


INTRODUCTION 

The  objective  of  this  study  was  to  calculate  the  smooth  wall,  thermo- 
chemical,  steady-state  ablation  response  of  selected  reentry  vehicle  nose  tip 
models  tested  as  part  of  the  SAMSO,  irospace  Thermal  Protection  Program.  Cal¬ 
culations  were  performed  for  POCO  graphite  and  R-6300  carbon-phenolic  models 
tested  in  the  Cornell  Wave  Superheater,  the  AFFDL  50  MW  RENT  facility,  and  the 
AEDC  Hypervelocity  Range  G  facility.  These  calculations  speciiically  exclude 
the  effects  of  the  rough  model  surface  on  heat  and  mass  transfer  rates  and  the 
effects  of  mechanical  ablation  on  the  model  response.  Thus,  these  predictions 
provide  an  accurate  "baseline"  for  comparison  to  experimentally  observed  abla¬ 
tion  response,  and  these  comparisons  provide  useful  input  to  mechanical  abla¬ 
tion  and  surface  roughness  effects  studies. 

The  test  conditions  and  model  contours  for  the  ablation  response  pre¬ 
dictions  were  taken  from  References  1  and  2  for  models  tested  in  the  CAL  WSH, 
from  Reference  3  for  models  tested  in  the  RENT  facility,  and  from  References 
4  and  5  for  models  tested  in  the  ballistic  range.  The  selection  of  the  particu¬ 
lar  models  for  analysis  was  guided  by  the  desire  to  consider  (1)  a  variety  of 
materials,  (2)  a  variety  of  test  facilities  and  conditions,  and  (3)  relatively 
high  quality  test  data,  particularly  contour  histories.  "Baseline"  smooth- 
wall  thermochemical  ablation  response  predictions  were  performed  for  the 
following  models: 


Model  No.  Test  Facility  Material 


G-9 

WSH 

POCO  graphite 

CP-16 

WSH 

R-6300  carbon 

phenolic 

12-2 

RENT 

POCO  graphite 

10  3A 

RENT 

R-6300  carbon 

phenolic 

123A 

RENT 

R-6300  carbon 

phenolic 

1750 

Ballistic  Range 

POCO  graphite 

1756 

Ballistic  Range 

POCO  graphite 

1763 

Ballistic  Range 

R-6300  carbon 

phenolic 

1861 

Ballistic  Range 

R-6300  carbon 

phenolic 
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In  addition  to  performing  "baseline"  thermochemical  ablation  calcula¬ 
tions,  a  number  of  parametric  calculations  for  one  of  the  models  tested 
in  the  wave  superheater  were  performed  in  order  to  quantitatively  assess 
the  effects  of  several  uncertainties  and  aasnmpt < nn«  on  i-So  oomnnfrod 
ablation  responses.  The  ur.certainties  assessed  in  this  fashion  were: 

o  uiicei Laintiea  associated  with  defining  the  model  contour  from  the 
reduced  movie  film  data  -  nose  radius  and  bi-conic  cone  angle 
uncertainties 

o  Transition  location  uncertainties 

o  The  effect  of  the  pressuiv  distribution  representation 

o  The  effect  of  unequal  spL-cies  diffusion  coefficients  for  carbon 

phenolic  ablation 

o  The  effect  of  the  "enthalpy  spike"  in  the  AFFDL  50  MW  RENT  facility. 

The  analytical  techniques  employed  in  these  computations  are  briefly 
described  in  Section  2.  The  baseline  model  response  predictions  are  compared 
to  measurements  in  Section  3,  and  appropriate  conclusions  are  drawn  in  Sec¬ 
tion  4, 


SKCTION  2 


ABLATION  AND  STRUCTURAI.  ANAr.VRTK  Cf^MPUTEP.  CODEE 

In  this  section,  the  computerised  ana]ytic<i]  tools  utilized  in  this 
program  are  described.  These  codes  evaluate  the  fluid  dynamics,  boundary 
layer  mass  and  heat  transfer,  surface  thermochemistry,  transient  in-depth  heat 
conduction,  and  thermal-structural  response  phenomena.  The  various  codes  are 
described  in  the  following  paragraphs. 

AEROTHERM  CHEMICAL  EQUILIBRIUM  (ACE)  CODE 

The  ACE  code  evaluates  the  equilibrium  or  kinetically  controlled  chemical 
composition  of  a  set  of  molecular  species  and  associated  thermodynamic  and  trans¬ 
port  properties  for  two  types  of  systems;  closed  systems  and  open  systems.  In 
a  closed  system  calculation,  the  relative  amounts  of  each  element  in  the  total 
system  plus  system  pressure  and  either  system  entropy,  temperature,  or  enthalpy 
are  specified.  For  an  open  system  calculation,  the  relative  amounts  of  the 
chemical  elements  depend  upon  various  boundary  layer  mass  transfer  formulations 
and  upon  the  specified  injection  rates  of  solid  material  decomposition  products 
and  surface  ablation  products.  Chemical  equilibrium  with  a  solid  surface 
species  provides  an  additional  constraint  in  open  system  computations.  Further 
descriptions  of  ACE  code  options  are  included  in  References  6,  7  and  8. 

INVISCID  FLOW  CODES,  RABBI  AND  SUPER 

Surface  pressure  distributions  and  bow  shock  shapes  are  obtained 
utilizing  two  computer  codes,  one  for  the  solution  of  the  subsonic  and  tran¬ 
sonic  nose  region,  (Realgas  Axisymmetric  Blunt  Body  at  Incidence,  Refs.  9  and 
10),  and  the  other  for  the  solution  of  the  supersonic  regions.  (SUPER,  Ref. 11) 
Solutions  in  the  subsonic/transonic  region  performed  in  this  study  were 
obtained  u.sing  a  zero  angle  of  attack  version  of  RABBI  code  described  in 
Reference  10.  The  steady  flow  field  was  obtained  as  the  time  asymptotic  limit 
of  an  unsteady  flow.  The  shock  was  treated  as  a  sharp  discontinuity.  The  code 
is  applicable  to  thermally  and  caloricalJy  perfect  gases  or  to  air  in  chemical 
equilibrium  (perfect  gas  equation  of  state) . 

The  supersonic  flow  field  program  (SUPER)  employs  the  method  of 
characteristics  as  described  in  Reference  11.  All  shocks  are  treated  as  sharp 
discontinuities;  sharp  expansion  corners  are  treated  as  concentrated  centered 
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expansions;  the  fluid  is  treated  as  a  thermally  and  calorically  perfect  gas. 

The  initial  data  for  this  code  is  obtained  from  the  RABBI  code,  described  in 

the  pr^  CCidii..^  ^ui.  a<.j  1  apli  . 

In  order  to  make  a  consistent  coupling  between  the  two  codes,  the 
solution  in  the  subsonic/transonic  region  is  done  for  an  ideal  gas  (here,  y=1.2 
was  used.) .  The  solutions  were  compared  with  the  solutions  for  air  in  equilib¬ 
rium,  and  it  was  found  that  there  was  no  difference  in  surface  pressure  distri¬ 
bution  and  negligible  difference  in  the  bow  shock  shape. 

BOUNDARY  LAYER  INTEGRAL  MATRIX  PROCEDURE  (BLIMP) 

The  BLIMP  code  computes  numerical  solutions  of  the  nonsimilar  multi- 
component  laminar  and/or  turbulent  boundary  layer  with  arbitrary  equilibrium 
or  nonequilibrium  chemically  reacting  systems  over  planar  or  axisymmetric  bodies. 
A  variety  of  surface  boundary  conditions,  including  coupling  with  steady  state 
or  transient  charring  ablation  energy  and  mass  balances,  are  available.  Multi- 
component  diffusion  and  thermal  diffusion  are  treated  through  the  use  of 
convenient  correlation  equations.  The  code  computes  its  own  boundary  layer 
edge  condition  from  the  stagnation  conditions,  the  pressure  distribution,  and, 
if  entropy  layer  is  being  considered,  shock  shape  input  data. 

The  numerical  solution  method,  termed  an  integral-matrix  procedure 
and  description  References  12,  13  and  14,  incorporates  sets  of  connected 
quadratics  between  boundary  layer  nodal  points  to  relate  enthalpy,  velocity, 
and  elemental  mass  fraction  variations  to  their  derivatives  with  respect  to 
the  boundary  layer  thicl<.ness  coodinate.  The  derivatives  of  these  quadratics 
are  made  continuous  at  the  nodal  points  thus  resulting  in  smooth  but  flexible 
profiles.  The  turbulent  flow  model  built  into  the  BLIMP  code  consists  of 
dividing  the  flow  into  a  wall  region  in  which  a  mixing  length  description 
of  turbulent  shear  is  used  and  a  wake  region  in  which  eddy  viscosity  is 
related  to  global  parameters  of  the  flow.  Further  discussion  of  the  BLIMP 
code  including  descriptions  of  the  turbulent  model  are  included  in  Reference 
15,  16  and  17. 

CHARRING  MATERIAL  ABLATION  (CMA)  PROGRAM 

The  CMA  program  computes  the  in-depth  transient  thermal  response  of  a 
one-space  dimension  charring,  ablating  material.  An  Arrhenius  rate  law  is 
utilized  to  obtain  material  decomposition  in-depth.  The  program  is  an 
implicit,  finite-difference  computational  procedure  for  heat  conductions  coupled 
to  ono  of  three  ablating  surface  boundary  conditions.  These  options  are: 
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•  A  film  coefficient  model  which  accounts  for  convection  - 
radiation  heating  with  coupled  mass  transfer,  including 
the  effects  of  unoaual  heat  and  mass  tranqfer.  iin»»giiai 
species  diffusion  rates,  and  surface  chemical  reactions 
in  equilibrium  or  kinetically  controlled. 

•  Specified  surface  temperature  and  recession  rate. 

•  Specified  radiation  view  factor  and  incident  radiation 
flux  for  a  stationary  surface. 

For  the  first  and  last  options  a  surface  energy  balance  is  utilized  in 
conjunction  with  the  in-depth  implicit  finite  difference  solution  to  define 
surface  temperature.  Surface  temperature  is  specified  in  the  second  option. 
Surface  thermochemistry  data  generated  by  the  ACE  code  is  utilized  by  the 
first  option.  Additional  descriptions  of  the  use  and  operation  of  the  CMA 
code  are  Included  in  References  18,  19  and  20. 

AXISYMMETRIC  TRANSIENT  HEATING  AND  MATERIAL  ABLATION  (ASTHMA)  PROGRAM 

The  ASTHMA  code  is  a  two  apace  dimension  transient  heat  conduction  code 
for  ablating  but  non-charring  materials.  It  accounts  for  two-dimensional 
surface  recession,  and  allows  for  anisotropic  in-depth  thermal  conduction. 

Surface  boundary  conditions  are  in  general  the  same  as  those  described  above 
for  the  CMA  code.  Input  to  the  ASTHMA  code  consists  of  geometry  specifications 
(in  terns  of  a  finite  difference  grid),  material  thermal  properties  tables,  and 
heated  surface  boundary  condition  information  as  a  function  of  lime,  output 
consists  of  two  dimensional  thermal  properties  ans  surface  recession.  Provision 
is  made  for  subsequent  input  of  thermal  data  into  various  structural  analyses 
codes.  Further  description  of  the  code  are  included  in  Reference  21  and 
addenda  thereto. 

DEFORMATION  FORMULATION  OF  ORTHOTROPIC  AXISYMMETRIC  SOLUTION  OF  INELASTIC 
SOLIDS  (DOASIS) 

The  DOASIS  finite  element  stress  analysis  code  determines  the  anisotropic, 
elastic-plastic-thermal  .solution  (i.e.  displacements,  stresses,  strains,  etc.)  of 
any  two  dimensional  structure  (axisymmetric,  plane  stress  or  plane  strain)  for 
any  prescribed  mechanical  and/or  thermal  load.  The  inelastic  response  is  based 
on  deformation  plasticity  theory,  in  conjunction  with  the  appropriate 
orthotropic  material  formulation.  The  finite  element  method  of  solution  is 
employed  which  subdivides  the  structure  into  a  number  of  "finite  elements"  with 
nodes.  The  solution  consists  of  finding  the  displacements  of  these  nodes  and 
then  the  resulting  strains  and  stresses  in  each  finite  element.  When  inelastic 
phenomenon  occurs,  the  nonlinear  solution  is  found  by  iteration.  Further  details 
of  the  theory  and  methodology  of  the  DOASIS  code  may  be  found  in  References 
22  and  23. 
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SECTION  3 


MATERIAL  RESPONSE  PREDICTIONS 

This  section  presents  experimental  data,  describes  how  the  baseline 
calculations  were  performed,  and  compares  predicted  to  measured  model  response. 
Calculations  for  the  ballistic  range  experiments  are  presented  in  Section  3.1 
and  are  followed,  in  Section  3.2  by  a  description  of  calculations  for  the  wave 
superheater  and  the  50  MW  RENT  facility. 

3.1  BALLISTIC  RANGE  EXPERIMENTS 

The  recent  improvements  in  the  instrumentation  and  data  acquisition 
capabililities  of  the  AEDC  1000  ft  Ballistic  Range  facility  have  enabled  ex¬ 
perimental  data  to  be  obtained  for  conditions  closely  approximating  those 
encountered  by  high  performance  reentry  vehicles.  These  data  should  be  useful 
not  only  in  experimentally  verifying  the  applicability  and  survival Llity  of 
new  heat  shield  materials,  but  also  in  the  evaluation  and  improvement  of  analyt¬ 
ical  techniques  for  the  prediction  of  material  response.  In  this  section,  de¬ 
tailed  predictions  of  the  thermal  response  of  four  model  shots  in  the  AEDC 
Range  G  facility  during  the  recently  completed  SAMSO/Aerospaoe  test  program^ 
arc  described.  Comparisons  between  predictions  and  test  measurements  provide 
a  basis  for  assessing  the  quality  of  the  data  and  the  current  analytical  tech¬ 
niques,  The  predictions  were  performed  utilizing  the  computer  codes  described 
above  in  Section  2.  The  test  conditions  of  the  range  shots  for  which  the  cal¬ 
culations  correspond  are  given  in  Section  3.1,1.  The  inviscid  flow  solution 
(shock  shape  and  pressure  distribution)  are  presented  in  Section  3.1.2.  Re¬ 
sults  of  boundary  layer  computations  are  discussed  in  Section  3.1.3.  The  sur¬ 
face  and  in-depth  transient  thermal  response  predictions  (including  recession) 
are  shown  in  Section  3.1.4.  And  the  evaluation  of  model  thermal  expansion  dur¬ 
ing  one  of  the  range  shots  is  described  in  Section  3.1,5. 

3.1.1  Test  Conditions 

Four  ballistic  range  shots  were  selected  for  analysis  in  this  study. 

The  tu'st  conditions,  model  geometries  and  materials  are  summarized  in  Table  I. 
Ml  four  shots  were  rformed  in  the  AEDC  1000  ft  ballistic  range  [^during  the 
SAHl'O/Ai-'rospac-i’  tost  ^■i(K)ruii.  Jtwo  shots  for  each  of  two  materials  were  chosen 
Lor  these  coacul.iti  one  to  demonstrate  the  effect  of  both  different  materials 
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SUMMARY  OF  BALLISTIC  RANGE  SHOTS 


SHOT  No. 

1750 

1756 

1763 

1861 

Material  (Layup  Angle) 

POCO 

POCO 

R-6300  (30°) 

R-6300(30°) 

Range  Temperature  (°F) 

75° 

75° 

75° 

75° 

Range  Static  Pressure  (atm) 

0.364 

0.458 

0.366 

0.236 

Blast  Tank  Static  Pressure (atm) 

0.013 

0.037 

0 . Oli' 

Launch  Velocity  (ft/sec) 

17000 

17900 

18580 

Time  to  Last  Station? 

(millisec) 

57.9 

62.1 

56.9 

52.2 

Initial  Geometry (Sphere-Biconic) 
Nose  Radius (inches) 

0.040 

0.040 

0.040 

0.045 

Fore  Cone  Half  Angle  (Deg.) 

45° 

45° 

40° 

37.5° 

Aft  Cone  Half  Angle  (Deg.) 

16° 

16° 

16° 

16° 

Ballistic  Coefficient (psf) 

58 

58 

58 

78 

Init.stag.pt. Pressure^  (atm) 

122.1 

139.2 

123.3 

85.8 

Init.  Stag.  Pt.  Cnthalp;^ 
(Btu/lbm) 

6372 

5773 

6400 

6896 

Stag.Pt.  Recession  (mils) 

25 

24 

44 

36 

(a)  Time  from  launch,  Station  42(940  ft.,  including  bJast  tank) 

(b)  Entrance  to  range  tank. 

(c)  JANAF  base  state,  298°K. 

(d)  Station  42. 


■  uui  varyinij  toHt  conditions.  In-f’iy!it  model  profile  traces  and  the  correspond- 
ui‘i  rcjosnidii  ineasuroiuents  wore  studied  to  define  those  shots  which  provided 
the  most  con;;  intent  set  of  response  data.  Shots  were  graded  not  only  on  the 
clarity  and  resolution  of  the  stagnation  point  recession  measurements,  but  also 
on  the  degree  of  symmetry  and  quality  of  shape  data  on  the  fore-cone  section 
cf  th^  ...odd.  i’leiiee  data  fui.  tile  siiots  predicted  shouia  provide  valuable  com¬ 
parisons  for  both  the  laminar  nose  region  and  the  turbulent  cone  region.  The 
tv.'c  matcri.alo  for  which  predictions  were  petfurineu  (FOCO  graphite  and  K-630U 
carbon  phenolic)  are  representative  of  currently  utilized  reentry  vehicle  ther¬ 
mal  protection  systems.  Figure  I  shows  the  overall  shape  of  the  two  POCO 
graphite  models. 

Model  velocity  histories  determined  from  test  data  for  each  shot  are 
shown  in  Figure  2.  The  velocity  decay  is  due  to  aerodynamic  drag  on  the  test 
model.  Based  upon  real  gas  computations  of  flow  past  a  normal  shock  wave 
(ACE  code) ,  the  stagnation  point  enthalpy  and  pressure  histories  were  obtained 
from  the  range  static  conditions  and  the  velocity  histories.  Stagnation  point 
pressure  and  enthalpy  histories  are  shown  in  Figures  3  and  4,  respectively, 
for  each  range  shot  considered.  These  histories  along  with  the  inviscid  flow 
solution  discussed  in  the  following  section  provide  the  necessary  input  to  the 
boundary  layer  and  thermal  structural  calculations. 

3.1.2  Inviscid  Flow  Solutions 

The  inviscid  flow  about  the  ballistic  range  model  is  one  of  the  primary 
inputs  to  the  evaluation  of  the  boundary  layer  flow  and  the  model  thermal- 
structural  response.  An  inviscid  flow  solution  for  a  particular  body  geometry 
consists  of  the  bow  shock  wave  geometry  and  the  associated  pressure  gradients. 

To  evaluate  these,  the  RABBI  and  SUPER  codes  discussed  in  Section  2  were  uti¬ 
lized  assuming  ideal  gas  flow  (y  =  1.2),  For  the  hypersonic  conditions  of 
interest  here  (Mach  numbers  between  11.0  and  16.5),  effects  of  varying  stagna¬ 
tion  pressure  and  Mach  number  on  shock  shape  and  pressure  distribution  are 
negligible.  This  was  verified  by  performing  calculations  at  various  free 
stream  conditions  for  the  45°  half-angle  sphere-cone  nose  geometries  used  in 
shots  1750  and  1756  and  shown  in  Figure  1. 

The  variation  of  model  geometry  due  to  ablation  during  flight  does  alter 
the  inviscid  flow.  To  account  for  this  variation,  in-flight  model  profile 
trace.s  wore  utilized  to  define  approximate  shapes  for  each  time  at  which  bound¬ 
ary  layer  solutions  were  obtained.  Because  of  the  extremely  short  exposure 
times  (-60  milliseconds),  recession  was  small  so  that,  although  ablation 
occurred,  the  models  remained  roughly  sphere-cones.  That  is,  ablation  resulted 
in  only  slight  increases  in  the  nose  radii  and  for  the  R-6300  models  slight  de¬ 
creases  in  the  cone  half  .xngle.  As  a  result,  sphere-cone  solutions  of  the 
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inviscid  flow  could  be  used  throughout.  This  provided  a  great  simplification 
since  inviscid  flow  solutions  are  identical  for  any  similar  geometrical  shape. 
Figure  S  shvjwa  Lhe  variation  ot  pressure  ratio  with  normalized  stream  dimension 
and  Figure  6  gives  the  shock  angle  as  a  function  of  normalized  radius  (from 
the  body  centeiline)  to  the  shock.  The  solution  for  the  45“^  half  angle  cone 
was  Used  for  all  POCO  graphite  model  solutions  while  the  solutions  between  35° 
and  40°  were  utilized  for  the  various  geometries  obtained  by  the  R-63t)0  models. 

3.1.3  Boundary  Layer  Solutions 

The  BLIMP  code  described  in  Section  2  was  utilized  to  compute  the  heat 
transfer  coefficient  distributions  around  the  ballistic-range  models  for  sev¬ 
eral  times  during  each  shot  analyzed.  The  solutions  accounted  for  the  varia¬ 
tions  in  free  stream  conditions,  model  geometry,  and  shock  shape  as  described 
above,  cross  plots  of  transfer  coefficient  versus  time  for  the  body  locations 
of  interest  provided  the  coefficient  histories  required  as  input  to  the  thermal 
response  codes . 

The  BLIMP  solutions  for  both  the  POCO  and  the  R-6300  models  were  per¬ 
formed  utilizing  a  surface  energy  balance  formulation  which  accounts  for  steady 
state  ablation  of  the  nose  tip  material  in  thermochemical  equilibrium  with  the 
gas  at  the  surface.  Also,  the  assumption  of  equal  diffusion  coefficients  for 
all  species  in  the  boundary  layer  wa.s  utilized.  This  assumption  is  further 
discussed  in  Section  3.2  belov/.  No  attempt  was  made  to  account  for  surface 
roughness  in  these  boundary  layer  calculations. 

The  effect  of  the  large  entropy  gradient  at  the  boundary  layer  edge 
due  to  the  shock  wave  curvature  was  quite  significant  in  these  soluti.ons.  This 
entropy  gradient  occurs  because  the  mass  at  different  locations  in  the  boundary 
layer  and  the  edge  gas  has  entered  the  flow  through  different  shock  angles.  In 
regions  where  the  shock  curvature  changes  rapidly,  this  entropy  layer  results 
in  an  edge  vorticity  much  different  than  zero.  The  effect  on  the  predicted 
heat  transfer  coefficient  distribution  is  shown  in  Figure  7.  Several  other 
solutions  for  the  same  geometry  are  shown  for  comparison.  As  would  be  expected, 
the  entropy  layer  solution  near  the  stagnation  point  is  near  the  solution  found 
using  only  a  normal  shock  entropy,  while  out  on  the  conic  region,  the  entropy 
layer  solution  approaches  that  obtained  using  sharp  cone  shock  angle  to  define 
the  edge  entropy.  Velocity  profiles  at  a  streamwise  location  of  0.008  feet 
for  the  three  distributions  are  shown  in  Figure  8.  In  the  figure,  edge  velocity 
is  used  to  normalize  the  local  velocity  in  the  boundary  layer.  The  effect  of 
the  vorticity  (i.e.,  the  entropy  layer)  is  to  produce  the  non-zero  slope  in 
profile  at  the  boundary  layer  edge. 
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The  minimum  point  in  the  coefficient  distribution  (at  a  streamwise  loca¬ 
tion  of  0.010  feet)  of  Figure  7  occurs  because  the  edge  vorticity  passes  through 
a  maximum  at  this  point.  The  velocity  gradients  within  the  boundary  layer  are 
least  severe  and  the  minimum  in  coefficient  results.  Some  question  exists  as 
to  the  reality  of  this  edge  vorticity  effect  on  the  coefficient  distributions. 
During  the  running  of  the  various  BLIMP  solutions  considerable  difficulty  was 
encountered  in  several  cases  at  this  location  on  the  body.  Convergence  was 
often  defeated  by  apparent  instabilities  in  the  numerical  technique  or  in  the 
turbulent  model  being  employed.  At  times,  more  than  one  solution  appeared 
mathematically  possible.  Once  past  this  region,  the  solutions  were  quickly 
obtained  and  appeared  correct.  As  seen  in  Figure  7,  downstream  of  the  minimum 
the  solution  approaches  that  for  the  sharp  cone  since  edge  vorticity  decreases 
to  zero  and  edge  entropy  approaches  the  sharp  cone  shock  angle  value. 

A  feeling  for  the  development  of  the  boundary  layer  can  be  obtained 
from  Figure  9.  Shown  are  five  velocity  profiles  for  various  locations  on  the 
body.  Local  velocity  in  this  figure  is  normalized  by  a  reference  velocity 
(the  velocity  which  is  calculated  for  an  isentropic  expansion  to  the  local 
pressure  reatio  for  a  normal  shock) .  Notice  that  due  to  entropy  layer  this 
reference  velocity  can  be  leas  than  the  computed  edge  velocity.  The  five  pro¬ 
files  shown  correspond  to  the  following  body  locationr?: 

•  A  -  Laminar  profile  near  sonic  point 

•  B  -  First  turbulent  profile 

•  C  -  Solution  upstream  of  vorticity  maximum  sh.owing  entropy  layer 

effects 

•  D  -  Solution  downstream  of  vorticity  maximum  showing  boundary  layer 

thinning  due  to  entropy  layer  effects 

•  E  -  Last  solution  approaching  sharp  cone  profile. 

The  heat  transfer  coefficient  distributions  for  each  BLIMP  solutior  per¬ 
formed  in  this  study  are  shown  in  Figure  10.  The  trends  discussed  aboVe  oocur 
throughout  although  the  locations  and  degree  of  severity  change  be  ■■.uk  -  m'  the 
varying  test  conditions  and  geometries.  For  the  two  carbon  phenolj.  s)  pre¬ 
dictions  were  carried  out  only  to  the  middle  of  the  fore  cone  where  i.^rit 

CMA  calculation  was  performed.  Examples  of  typical  boundary  layer  co'''position 
profiles  are  shown  in  Figure  11.  One  is  for  a  POCO  graphite  prediotioti  and 
one  is  for  an  R-6300  carbon  phenolic  prediction. 
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3.1.4  Comparison  of  Measured  and  Predicted  Ablation  Responses 

The  thermal  and  ablation  response  predictions  ror  the  two  ruCO  yraptiite 
shots  were  performed  by  the  ASTHMA  code.  The  CMA  code  was  utilized  to  compute 
the  thermal  and  ablation  response  for  two  locations  on  each  of  the  two  R-6300 
carbovi  phenolic  models  of  interest.  The  ACE  code  was  used  to  generate  the 
equilibrium  surface  thermochemistry  data  required  as  input  to  both  the  CMA  and 
ASTHMA  codes.  Plots  of  some  of  these  data  are  given  in  Figure  12.  Shown  for 
the  POCO  graphite  predictions  are  the  normalized  ablation  rate  versus . surface 
temperature  for  various  values  of  pressure.  For  the  R-6300  carbon  phenolic, 
normalized  ablation  rate  versus  surface  temperature  for  various  noiii.alized 
pyrolysis  gas  rates  at  a  particular  pressure  are  shown. 

Material  thermal  conductivity  variations  with  temperature  are  given  in 
Figure  13.  The  values  used  in  the  R-6300  predictions  were  those  of  a  very 
similar  material  MX-4926  carbon  phenolic  at  a  0°  layup  angle.  Extrapolation 
to  temperatures  above  8000°R  was  required. 

Various  results  of  the  transient  thermal  response  predictions  are  given 
in  Figures  14  through  16.  Figure  14  shows  the  variations  of  surface  tempera¬ 
ture  at  various  locations  on  the  respective  models.  Also  shown  as  symbols  are 
various  steady  state  surface  temperatures  computed  by  the  BLIMP /»ode.  This 
comparison  verifies  that  essentially  steady  state  ablation  occurs  throughout 
the  ballistic  range  shots. 

Figure  15  gives  several  typical  predicted  in-depth  temperature  and,  for 
the  R-6300  carbon  phenolic,  density  profiles  at  the  last  measurement  station. 

The  figure  indicates  that  the  back  wall  temperatures  in  the  POCO  graphite 
models  are  high  enough  to  cause  the  epoxy  bonds  and  Sootchply  backup  to  begin 
decomposing  (-lOOD^R) .  Heat  soak  did  not  reach  the  back  surface  of  the  R-6300 
mode Is ,  however . 

The  predictions  of  nose  region  recession  for  the  four  shots  are  shown 
in  Figure  16.  Comparison  is  made  between  the  final  measured  geometries  and 
the  predicted  surface  locations  data  points.  For  shot  1750,  an  additional 
prediction  which  accounts  for  thermal  expansion  is  shown.  This  prediction 
will  be  discussed  in  the  following  section.  In  general,  less  ablation  was 
predicted  than  was  indicated  by  the  in-flight  model  photographs.  At  the  stag¬ 
nation  point  for  all  cases,  predicted  recession  was  about  1/3  of  measured, 
while  on  the  cone  region,  predicted  was  only  about  1/2  of  measured. 

3.1.5  Thermal  Expansion  Calculations 

Two  elastic  thermal  expansion  calculations  were  made  utilizing  the  DOASIS 
code  described  in  Section  2.  The  ASTHMA  prediction  of  shot  1750  prov.  -ad  the 
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necessary  thermai  and  mechanical  loading  plus  the  geometrical  input  data.  The 
finite  element  network  and  temoerature  distributions  employed  are  shown  in  Fig¬ 
ures  17a  and  b  respectively.  Thermal-structural  properties  of  the  POCO  graphite 
were  obtained  from  Reference  24.  Predictioiis  were  made  for  two  back  wall  sup¬ 
port  conditions.  These  were: 

•  Nose  tip  supported  by  Scotchply  plug  on  the  back  wall  (Surface  A 
in  Figure  1) 

•  Nose  tip  supported  by  beryllium-copper  afterbody  at  the  end  of  the 
aft-cone  [Surface  B  in  Figure  1) 

The  second  condition  more  accurately  represents  the  actual  support  condition 
for  shot  1750.  For  comparison  purposes,  both  cases  were  run.  For  the  support 
at  surface  A,  the  thermal  expansion  at  the  stagnation  point  was  2.5  mils  and 
less  than  1  mil  at  the  bi-conic  intersection.  For  the  support  at  surface  B, 
the  thermal  expansion  was  6.5  mils  at  the  stagnation  point  and  4.0  mils  at  the 
bi-conic  intersection.  The  differences  between  these  two  surfade  expansions 
comes  from  the  differential  thermal  expansion  between  support  A  and  B.  Figure 
16a  shows  the  predicted  final  contour  corrected  for  the  latter  thermal  expansion. 

Sufficient  high  temperature  material  properties  data  was  not  available 
to  adequately  perform  a  reliable  plasticity  analysis,  however,  it  should  be 
pointed  out  that  inclusion  of  plastic  flow  at  the  outer  surface  (T  >  5000°F) 
would  decrease  the  predicted  thermal  expansion.  Consequently,  the  thermal  ex¬ 
pansion  shown  in  Figure  16a  should  be  considered  as  an  upper  limit. 

3.2  WAVE  SUPERHEATER  AND  50  MW  RENT  FACILITY  TESTS 

Calculations  presented  above  for  ballistic  range  tests  required  consid¬ 
eration  of  the  changing  free  stream  conditions  and  to  some  extent,  changing 
model  contour  with  time.  The  wave  superheater  and  RENT  facility  tests  are 
characterized  by  only  minor  changes  in  test  conditions  with  time  and  the  par¬ 
ticular  test  models  chosen  for  analysis  were  selected  on  the  basis  of  having 
established  a  relatively  stable  shape.  Because  of  the  high  convective  heating 
rate  conditons,  it  is  reasonable  to  assume  that  steady  state  energy  and  mass 
balances  at  the  surface  are  representative  of  reality  for  test  times  of  approx¬ 
imately  1  second  and  longer.  Hence,  for  each  model  analyzed  an  ablation  predic¬ 
tion  was  performed  for  one  point  in  time.  The  test  conditions  and  model  con¬ 
tours  are  described  first,  in  Section  3.2.1,  and  are  followed  by  a  description 
of  inviscid  flow  representation  and  boundary  layer  solutions  in  Sections  3.2.2 
and  3.2.3  respectively.  Predicted  ablation  response  and  comparison  to  data 
are  given  in  Section  3.2.4. 
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Test  Conditions  and  Model  Contours 


For  the  models  tested  in  the  WSH,  the  stagnation  point  impact  pressure 
and  total  enthalpy  (Hq)  values  used  in  the  calculations  were  taken  di¬ 
rectly  rrom  Reterence  1.  The  enthalpy  values  reported  in  Reference  1  were 
decreased  by  129.  Btu/lbm  to  be  consistent  with  the  JANAF  enthalpy  base  state 
employed  nere. 

For  the  models  tested  in  the  RENT  facility,  the  stagnation  point  pres¬ 
sures  used  in  the  analyses  accounted  for  the  flight  impact  pressure  decay  along 
the  stream  centerline.  At  a  given  axial  location  the  impact  pressure  was  ob¬ 
tained  from  the  measured  chamber  pressure  (Ref.  3)  and  the  ratio  of  impact  to 
chamber  pressure  shown  in  Figure  18  taken  from  the  measurements  reported  in 
Reference  25.  The  bulk  stream  enthalpy  values  for  the  models  tested  in  the 
RENT  facility  are  reported  in  Reference  3.  However,  calorimetric  heat  flux 
measurements  reported  in  References  25  and  26  indicate  that  there  is  a  sig¬ 
nificant  centerline  "heat  flux  spike"  in  the  RENT  test  jet  at  high  pressure 
conditions.  Centerline  enthalpy  measurements  have  not  yet  been  made  for  the 
RENT  test  jet  although  heat  flux  measurements  and  heat  transfer  calculations 
in  Reft)rence  25  suggest  that  the  high  centerline  heat  flux  results  because 
the  centerline  enthalpy  is  on  the  order  of  5000  Btu/lbm  for  the  2400  Btu/lbm 
nominal  bulk  enthalpy  high  pressure  test  conditions.  Also,  the  studies  of 
Reference  27  demonstrate  that  the  model  response  is  essentially  unaffected 
by  the  enthalpy  of  the  fluid  away  from  the  vicinity  of  the  centerline  provided 
the  ablation  model  is  aligned  with  the  enthalpy  spike.  It  stands  to  reason 
that  a  misalignment  would  result  in  an  asymmetric  ablation  contour  and  this 
was  not  observed  for  the  models  that  were  analyzed.  Thus,  ablation  response 
calculations  for  models  tested  in  the  RENT  facility  utilized  an  enthlapy  of 
5000  Btu/lbm.  To  assess  the  affect  of  this  enthalpy  uncertainty,  an  additional 
calculation  was  carried  out  using  the  reported  bulk  enthalpy,  and  this  is  dis¬ 
cussed  in  Section  3.2.  The  test  conditions  considered  for  each  model  analyzed 
are  listed  in  Table  II. 

The  model  contour  histories  obtained  from  the  movie  film  data  and  the 
contours  analyzed  are  shown  in  Figures  19a  through  e.  The  film  data  was  ob¬ 
tained  from  References  1  and  3  and  the  model  contours  were  selected  to  repre¬ 
sent  an  average  effective  axisymirietric  steady  state  shape.  Figures  20a  through 
e  present  measured  surface  recession  as  a  function  of  time  for  all  models. 

Also  shown  in  these  figures  are  the  effective  cone  half  angle  and  nose  radii 
as  deduced  from  the  film  data  for  all  models  except  model  12-2  which  did  not 
assume  a  bi-conic  shape.  For  models  tested  in  the  wave  superheater  the  mea¬ 
sured  stagnation  point  surface  temperature  history  is  also  shown.  By  way  of 
example,  the  method  employed  to  select  the  time  at  which  analyses  were 
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performed  may  be  illustrated  by  referring  to  Figure  20a  for  model  G-9  tested 
in  the  wave  superheater.  For  this  particular  case,  the  contour  analyzed  was 
selected  to  have  a  64°  cone  half  angle  and  a  noae  i.<»uiuB  of  0.05  inch.  xnir. 
contour  corresponds  approxitnately  to  a  time  of  2.3  sec  which  is  after  most 
shape  change  has  occurrarf,  the  surface  temperature  has  leveled  out,  and  the 
recession  rate  is  relatively  constant. 

As  indicated  in  Figure  20a,  there  is  considerable  uncertainty  associated 
with  estimating  an  effective  cone  half  angle  and  nose  radius  from  the  film  data. 
WSH  tested  models  ablated  to  somewhat  asymmetric  shapes  as  indicated  in  Figures 
19a  and  b  while  models  tested  in  the  RENT  facility  (Figures  19c  through  19e) 
displayed  only  minor  asymmetries.  The  effects  on  the  ablation  predictions  of 
uncertainties  associated  with  defining  an  effective  nose  contour  from  the  movie 
film  data  were  assessed  by  carrying  out  additional  solutions  with  braeJeeting 
values  of  nose  radius  and  cone  half  angle.  It  was  found  that  the  effect  of 
these  uncertainties  on  predicted  recession  rate  was  from  10  to  15  percent. 
Detailed  results  from  these  solutions  are  given  in  Section  3.2. 

3.2.2  Inviscid  Flow  Representation 

Boundary  layer  solutions  require  specification  of  (1)  the  model  surface 
pressure  distribution,  vncl  (2)  if  shock  curvature  is  great,  the  shock  shape 
must  also  be  specified  in  order  to  properly  account  for  vortioity  at  the  bound¬ 
ary  layer  edge.  A  single  solution  was  performed  accounting  for  vorticity  and 
it  was  concluded  that  entropy  layer  effects  are  negligible  for  the  test  con¬ 
ditions  and  model  shapes  being  considered  here  (see  Table  II) .  Because  a  de¬ 
tailed  shock  shape  specification  was  not  required  it  was  not  necessary  to  per¬ 
form  numerical  inviscid  flow  field  solutions  as  was  required  for  the  ballistic 
range  tests  described  in  Section  3.1. 

Pressure  distributions  for  all  models  were  computed  employing  a  corrected 
form  of  the  pressure  distribution  correlation  suggested  in  Reference  28.  The 
correlation  suggested  in  Reference  28  gives  unsatisfactory  results  when  the 
free  stream  pressure  (nozzle  exit  pressure)  was  significant  relative  to  the 
model  stagnation  pressure.  The  pressure  distribution  correlation  equation 
employed  is  given  here. 
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is  the  maximum  of  either  R.,  or  R* 
max  T4 

*  denotes  sonic  condition 

Pi.  “  stagnation  point  impact  pressure 

P^  =  free  stream  static  pressure 

The  sonic  point  is  located  by  assuming  a  Newtonian  pressure  ratio  there  which 
yields: 
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Tho  above  pressure  distribution  equation  is  used  for  the  stagnation  point  to 
the  sonic  point.  A  Newtonian  P  given  ny 


P  =  P  +  (1  -  P  )cos^8 

or}  '  <jr)' 


is  used  downstream  of  the  sonic  point  and  up  until 


ds 


Newtonian 


ds 


Prandtl-Meyer 


and  P  is  obtained  from  the  Prandtl-Meyer  expansion  expression  downstream  of 
this  "match  point." 

Pressure  distributions  for  all  models  are  shown  in  Figures  21a  through 
e.  In  addition  to  the  pressure  distribution  for  the  baseline  contour  for  model 
G-9,  Figure  2ia  shows  pressure  distributions  employed  for  investigating  effects 
of  uncertainties  in  measurement  of  cone  angle  and  nose  radius  as  well.  Also 
shown  for  comparison  purposes  is  the  Newtonian  pressure  distribution  for  this 
model.  Boundary  layer  solutions  including  steady  state  surface  recession  rate 
calculations  were  performed  for  all  pressure  distributions  shown  in  Figures  21a 
through  e. 


3.2.3  Boundary  Layer  Solutions 

All  boundary  layer  solutions  were  performed  with  the  BLIMP  program.  With 
the  exception  of  one  solution  performed  for  model  G-9  accounting  for  entropy 
layer  effects,  all  solutions  were  performed  assuming  the  boundary  layer  edge 
thermodynamic  state  was  represented  by  an  equilibrium  isentropio  expansion 
from  stagnation  (normal  shock)  conditions.  For  all  models  which  assumed  a 
conical  nose  shape  (all  but  model  12-2  which  blunted)  boundary  layer  transi¬ 
tion  was  assumed  to  occur  just  past  the  sphere  cone  junction.  For  model  12-2 
(POCO  graphite  tested  in  the  RENT  facility)  the  boundary  layer  was  assumed 
laminar  over  the  entire  nose.  A  brief  study  was  also  conducted  to  assess  the 
effect  of  transition  location  upon  predicted  recession  for  the  models  which 
assumed  a  conical  shape.  Varying  the  transition  location  from  the  stagnation 
point  to  just  downstream  of  the  sphere  cone  junction  had  negligible  effect 
upon  predicted  recession  rate  in  the  mid-cone  region. 

All  solutions  but  one  employed  the  assumption  of  equal  molecular  diffu¬ 
sion  coefficients  and  the  neglect  of  thermal  diffusion.  One  solution  was  per¬ 
formed  with  best  estimate  molecular  diffusion  coefficients  for  all  species  and 
included  thermal  diffusion  as  well.  This  solution  was  performed  for  an  R-6300 
carbon  phenolic  model  (CP-16)  and  revealed  only  a  3  percent  effect  on  predicted 
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recession  rate  when  compared  to  the  equal  diffusion  calculation  for  the  same 
model . 

Steady  state  energy  and  mass  balances  were  assumed  for  all  models.  The 
steady  state  mass  balance  implies  that  the  relative  quantity  of  chemical  ele¬ 
ments  being  ablated  from  the  surface  are  in  the  same  ratio  as  they  are  in  the 
virgin  material.  The  steady  state  energy  balance  assumption  implies  that  the 
rate  of  heat  conduction  into  the  material  from  the  ablating  surface  is  equal 
to  the  product  of  mass  loss  rate  and  enthalpy  increase  of  the  ablation  material 
from  its  initial  temperature  to  the  surface  temperature. 

conduction  ~  *^initial^ ablation 

material 


3.2.4  Ablation  Response 


This  section  presents  results  directed  at  cissessing  the  effect  of  sev¬ 
eral  uncertainties  upon  predicted  ablation  rcsponue  and  compares  measured  to 
predicted  ablation  material  performance. 


3.2.4. 1  Effect  of  Uncertainties 

Prior  to  making  predictions  for  comparing  to  measured  data,  an.  assess¬ 
ment  was  made  of  the  effect  of  various  uncertainties  upon  predicted  ablation 
rate.  These  uncertainties  fall  Into  two  categories! 

1.  Experimental  Data  Uncertainties 

o  Inability  to  precisely  measure  nose  radius  from  film  data 

o  Inability  to  precisely  determine  the  body  cone  angle  from 

film  data 

o  Uncertainty  in  stream  total  enthalpy  in  RENT  facility 

2.  Uncertainties  in  Mathematically  Modeling  the  Physics 

o  Location  of  boundary  layer  transition 

o  Uncertainties  in  assessing  the  nose  pressure  distribution 

o  Are  entropy  layer  effects  important? 

o  Is  it  necessary  to  treat  unequal  molecular  diffusion 
coefficient  effects? 

With  the  exception  of  the  total  enthalpy  uncertainty  (2530  Btu/lb  bulk 
enthalpy  vs.  5000  Btu/lb  on  the  centerline  inferred  from  heat  flux  measurements) 
in  the  RENT  facility,  uncertainties  investigated  generally  resulted  in  less 
than  a  15  percent  change  in  surface  recession  rate  on  the  conical  (turbulent) 
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portion  of  the  test  model.  Utilization  of  the  bulk  enthalpy  rather  than  the 
estimated  centerline  value  in  the  RENT  facility  resulted  in  a  35  percent  de- 
crcazG  in  preulwLeu  ablation  rate. 

The  effect  of  many  of  the  investigated  uncertainties  upon  predicted 
ablation  rate  is  shown  in  Figures  22a  through  e.  A  tabular  summary  of  the 
effect  of  the  investigated  uncertainties  is  given  in  Table  III. 

3. 2. 4. 2  Comparison  of  Predictions  with  Data 

The  predicted  and  measured  axial  surface  recession  rates  are  shown 
in  Figures  23a  through  e.  The  axial  recession  rate  is  shown  in  the  figures 
for  comparison.  The  experimental  recession  rates  are  all  shown  as  horizontal 
lines  since,  at  the  times  analyzed,  relatively  stable  experimental  shapes  were 
observed.  Predicted  and  measured  stagnation  point  surface  teinperatures  are 
uIlo  shown  for  models  tested  in  the  wave  superheater. 

The  predicted  stagnation  point  surface  temperature  for  the  POCO  graphite 
model  in  the  wave  superheater  (Model  G-9)  is  about  350°R  above  the  measured 
value,  whereas  the  corresponding  prediction  for  carbon  phenolic  is  1300  R  above 
the  measurement,  It  is  believed  that  the  large  discrepancy  in  predicted  sur¬ 
face  temperature  for  the  carbon  phenolic  model  is  a  result  of  mechanical  erosion. 
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SECTION  4 


CONCLUSIONS 

The  study  objective  was  to  perform  baseline  calculations  of  ablative 
model  performance  considering  only  thermochemical  ablation  and  representing 
boundary  layer  heat-and-masa  transport  phenomena  by  smooth-wall  relations. 

The  principal  conclusion  to  be  reached  from  this  st.idy  relates  to  the 
discrepancy  between  these  baseline  predictions  and  observed  model  performance. 
Some  accessory  conclusions  are  also  reached  which  relate  to  the  effect  of 
various  uncertainties  upon  predicted  ablation  performance. 

Table  IV  summarizes  results  of  the  study  in  the  form  of  predicted-to- 
measured  recession  rate  ratio  for  the  model  stagnation  points  and  mid-cone 
sections.  This  ratio  is  plotted  as  a  function  of  local  surface  pressure  in 
Figure  -24.  With  the  exception  of  three  ballistic  range  stagnation  point  pre¬ 
dictions  all  data  correlate  well  with  local  pressure.  It  is  noted  that 
inclusion  of  thermal  expansion  effects  would  tend  to  bring  these  points  in  line 
with  the  correlation;  however,  some  of  the  other  ballistic  range  predictions 
would  not  correlate  as  well  if  thermal  expansion  effects  were  included.  They 
are  not  included  in  the  correlation  because  significant  uncertainty  exists  in 
the  magnitude  of  thermal  expansion  as  a  result  of  ill-defined  high  temperature 
property  data. 

It  is  concluded  that  surface  roughness  effects  and  mechanical  ablation 
effects  become  increasingly  important  as  pressure  increases.  The  choice  of 
pressure  as  the  correlation  parameter  here  is  primarily  in  the  interest  of 
expedience.  A  detailed  effort  to  correlate  the  data  in  terms  of  more 
fundamental  parameters  was  not  within  the  study  scope.  Other  parameters  which 
should  be  considered  include  surface  temperature,  rough  wall  shear  and  perhaps 
pressure  gradient. 

The  effects  of  various  experimental  uncertainties  upon  predicted  smooth 
wall  thermochemical  ablation  are  summarized  in  Table  III.  It  is  concluded  that 
uncertainties  in  transition  location  and  nose  radius  have  a  negligible  effect 
upon  predicted  turbulent  cone  ablation  rate.  Similarly  the  effects  of  unequal 
species  diffusion  coefficients  and  shock  layer  vorticity  have  negligible 
effects  for  models  tested  in  the  wave  superheater  and  RENT  facilities.  A  10  to 
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15  percent  error  in  predicted  cone  ablation  rate  is  associated  with  nominaJ 
uncertainties  in  cone  half  anqle  and  with  utilization  of  an  approximate 
(Newtonian)  pressure  distribution.  Vorticity  effects  are  very  important  over 
the  entire  conic  surface  for  models  tested  in  the  ballistic  range.  Because 
fundamental  data  is  lacking  upon  which  to  base  a  turbulent  prediction  with 
strong  edge  velocity  gradients  (normal  to  the  boundary  layer  edge)  these 
predictions  should  be  viewed  to  have  a  10  to  30  percent  uncertainly. 
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